Selective genotyping strategies for a sib test scheme of a broiler breeder program.

BACKGROUND: In broiler breeding, genotype-by-environment interaction is known to result in a genetic correlation between body weight measured in bio-secure and commercial environments that is substantially less than 1. Thus, measuring body weights on sibs of selection candidates in a commercial environment and genotyping them could increase genetic progress. Using real data, the aim of this study was to evaluate which genotyping strategy and which proportion of sibs placed in the commercial environment should be genotyped to optimize a sib-testing breeding program in broilers. Phenotypic body weight and genomic information were collected on all sibs raised in a commercial environment, which allowed to retrospectively analyze different sampling strategies and genotyping proportions. RESULTS: Accuracies of genomic estimated breeding values (GEBV) obtained with the different genotyping strategies were assessed by computing their correlation with GEBV obtained when all sibs in the commercial environment were genotyped. Results showed that, compared to random sampling (RND), genotyping sibs with extreme phenotypes (EXT) resulted in higher GEBV accuracy across all genotyping proportions, especially for genotyping proportions of 12.5% or 25%, which resulted in correlations of 0.91 vs 0.88 for 12.5% and 0.94 vs 0.91 for 25% genotyped. Including pedigree on birds with phenotype in the commercial environment that were not genotyped increased accuracy at lower genotyping proportions, especially for the RND strategy (correlations of 0.88 vs 0.65 at 12.5% and 0.91 vs 0.80 at 25%), and a smaller but still substantial increase in accuracy for the EXT strategy (0.91 vs 0.79 for 12.5% and 0.94 vs 0.88 for 25% genotyped). Dispersion bias was virtually absent for RND if 25% or more birds were genotyped. However, GEBV were considerably inflated for EXT, especially when the proportion genotyped was low, which was further exacerbated if the pedigree of non-genotyped sibs was excluded. CONCLUSIONS: When less than 75% of all animals placed in a commercial environment are genotyped, it is recommended to use the EXT strategy, because it yields the highest accuracy. However, caution should be taken when interpreting the resulting GEBV because they will be over-dispersed. When 75% or more of the animals are genotyped, random sampling is recommended because it yields virtually no bias of GEBV and results in similar accuracies as the EXT strategy.

Economic values of reproductive, growth, feed efficiency and carcass traits in Nellore cattle.

A bioeconomic model was developed to calculate the economic value (ev) of reproductive and growth performance, feed efficiency and carcass traits of a seedstock Nellore herd. Data from a full-cycle cattle operation (1,436 dams) located in the Brazilian Cerrado were assessed. The ev was calculated by the difference in profit before and after one-unit improvement in the trait, with others remaining unchanged. The ev was standardized by the phenotypic standard deviation of each trait. Preweaning average daily gain (ADG) was the most economically important trait evaluated (R$ 58.04/animal/year), followed by age at first calving (R$ 44.35), postweaning ADG (R$ 31.43), weight at 450 days (R$ 25.36), accumulated productivity (R$ 21.43), ribeye area (R$ 21.35), calving interval (R$ 19.97), feed efficiency (R$ 15.24), carcass dressing per cent (R$ 8.27), weight at 120 days (R$ 6.22), weight at 365 days (R$ 6.06), weight at weaning (210 days, R$ 5.82), stayability (R$ 5.70) and the probability of early calving (R$ 0.32). The effects of all traits on profits are evidence that their selection may result in the economic and genetic progress of the herd if there is genetic variability.

Weighted single-step genome-wide association study and pathway analyses for feed efficiency traits in Nellore cattle.

The aim was to conduct a weighted single-step genome-wide association study to detect genomic regions and putative candidate genes related to residual feed intake, dry matter intake, feed efficiency (FE), feed conversion ratio, residual body weight gain, residual intake and weight gain in Nellore cattle. Several protein-coding genes were identified within the genomic regions that explain more than 0.5% of the additive genetic variance for these traits. These genes were associated with insulin, leptin, glucose, protein and lipid metabolisms; energy balance; heat and oxidative stress; bile secretion; satiety; feed behaviour; salivation; digestion; and nutrient absorption. Enrichment analysis revealed functional pathways (p-value < .05) such as neuropeptide signalling (GO:0007218), negative regulation of canonical Wingless/Int-1 (Wnt) signalling (GO:0090090), bitter taste receptor activity (GO:0033038), neuropeptide hormone activity (GO:0005184), bile secretion (bta04976), taste transduction (bta0742) and glucagon signalling pathway (bta04922). The identification of these genes, pathways and their respective functions should contribute to a better understanding of the genetic and physiological mechanisms regulating Nellore FE-related traits.

Would large dataset sample size unveil the potential of deep neural networks for improved genome-enabled prediction of complex traits? The case for body weight in broilers.

BACKGROUND: Deep neural networks (DNN) are a particular case of artificial neural networks (ANN) composed by multiple hidden layers, and have recently gained attention in genome-enabled prediction of complex traits. Yet, few studies in genome-enabled prediction have assessed the performance of DNN compared to traditional regression models. Strikingly, no clear superiority of DNN has been reported so far, and results seem highly dependent on the species and traits of application. Nevertheless, the relatively small datasets used in previous studies, most with fewer than 5000 observations may have precluded the full potential of DNN. Therefore, the objective of this study was to investigate the impact of the dataset sample size on the performance of DNN compared to Bayesian regression models for genome-enable prediction of body weight in broilers by sub-sampling 63,526 observations of the training set. RESULTS: Predictive performance of DNN improved as sample size increased, reaching a plateau at about 0.32 of prediction correlation when 60% of the entire training set size was used (i.e., 39,510 observations). Interestingly, DNN showed superior prediction correlation using up to 3% of training set, but poorer prediction correlation after that compared to Bayesian Ridge Regression (BRR) and Bayes Cπ. Regardless of the amount of data used to train the predictive machines, DNN displayed the lowest mean square error of prediction compared to all other approaches. The predictive bias was lower for DNN compared to Bayesian models, across all dataset sizes, with estimates close to one with larger sample sizes. CONCLUSIONS: DNN had worse prediction correlation compared to BRR and Bayes Cπ, but improved mean square error of prediction and bias relative to both Bayesian models for genome-enabled prediction of body weight in broilers. Such findings, highlights advantages and disadvantages between predictive approaches depending on the criterion used for comparison. Furthermore, the inclusion of more data per se is not a guarantee for the DNN to outperform the Bayesian regression methods commonly used for genome-enabled prediction. Nonetheless, further analysis is necessary to detect scenarios where DNN can clearly outperform Bayesian benchmark models.

Genome-wide scan highlights the role of candidate genes on phenotypic plasticity for age at first calving in Nellore heifers.

Age at first calving (AFC) plays an important role in the economic efficiency of beef cattle production. This trait can be affected by a combination of genetic and environmental factors, leading to physiological changes in response to heifers' adaptation to a wide range of environments. Genome-wide association studies through the reaction norm model were carried out to identify genomic regions associated with AFC in Nellore heifers, raised under different environmental conditions (EC). The SNP effects for AFC were estimated in three EC levels (Low, Medium, and High, corresponding to average contemporary group effects on yearling body weight equal to 159.40, 228.6 and 297.6 kg, respectively), which unraveled shared and unique genomic regions for AFC in Low, Medium, and High EC levels, that varied according to the genetic correlation between AFC in different EC levels. The significant genomic regions harbored key genes that might play an important biological role in controlling hormone signaling and metabolism. Shared genomic regions among EC levels were identified on BTA 2 and 14, harboring candidate genes associated with energy metabolism (IGFBP2, IGFBP5, SHOX, SMARCAL1, LYN, RPS20, MOS, PLAG1, CHCD7, and SDR16C6). Gene set enrichment analyses identified important biological functions related to growth, hormone levels affecting female fertility, physiological processes involved in female pregnancy, gamete generation, ovulation cycle, and age at puberty. The genomic regions highlighted differences in the physiological processes linked to AFC in different EC levels and metabolic processes that support complex interactions between the gonadotropic axes and sexual precocity in Nellore heifers.

Investigating causal biological relationships between reproductive performance traits in high-performing gilts and sows1.

Efficient management of swine production systems requires understanding of complex reproductive physiological mechanisms. Our objective in this study was to investigate potential causal biological relationships between reproductive performance traits in high-producing gilts and sows. Data originated from a nutrition experiment and consisted of 200 sows and 440 gilts arranged in body weight blocks and randomly assigned to dietary treatments during late gestation at a commercial swine farm. Reproductive performance traits consisted of weight gain during late gestation, total number born and number born alive in a litter, born alive average birth weight, wean-to-estrous interval, and total litter size born in the subsequent farrowing. Structural equation models combined with the inductive causation algorithm, both adapted to a hierarchical Bayesian framework, were employed to search for, estimate, and infer upon causal links between the traits within each parity group. Results indicated potentially distinct reproductive networks for gilts and for sows. Sows showed sparse connectivity between reproductive traits, whereas the network learned for gilts was densely interconnected, suggesting closely linked physiological mechanisms in younger females, with a potential for ripple effects throughout their productive lifecycle in response to early implementation of tailored managerial interventions. Cross-validation analyses indicated substantial network stability both for the general structure and for individual links, though results about directionality of such links were unstable in this study and will need further investigation. An assessment of relative statistical power in sows and gilts indicated that the observed network discrepancies may be partially explained on a biological basis. In summary, our results suggest distinctly heterogeneous mechanistic networks of reproductive physiology for gilts and sows, consistent with physiological differences between the groups. These findings have potential practical implications for integrated understanding and differential management of gilts and sows to enhance efficiency of swine production systems.

Spatialization of climate, physical and socioeconomic factors that affect the dairy goat production in Brazil and their impact on animal breeding decisions / Espacialização de fatores climáticos, físicos e socioeconômicos que afetam a produção de caprinos leiteiros no Brasil e seu impacto sobre as decisões em melhoramento animal

Brazil has high climate, soil and environmental diversity, as well as distinct socioeconomic and political realities, what results in differences among the political administrative regions of the country. The objective of this study was to determine spatial distribution of the physical, climatic and socioeconomic aspects that best characterize the production of dairy goats in Brazil. Production indices of milk per goat, goat production, milk production, as well as temperature range, mean temperature, precipitation, normalized difference vegetation index, relative humidity, altitude, agricultural farms; farms with native pasture, farms with good quality pasture, farms with water resources, farms that receive technical guidance, family farming properties, non-familiar farms and the human development index were evaluated. The multivariate analyses were carried out to spatialize climatic, physical and socioeconomic variables and so differenciate the Brazilian States and Regions. The highest yields of milk and goat production were observed in the Northeast. The Southeast Region had the second highest production of milk, followed by the South, Midwest and North. Multivariate analysis revealed distinctions between clusters of political-administrative regions of Brazil. The climatic variables were most important to discriminate between regions of Brazil. Therefore, it is necessary to implement animal breeding programs to meet the needs of each region.

O Brasil possui diversidade edafoclimática e realidades socioeconômicas e políticas distintas. Isto contribui para diferenciar as regiões político administrativas do país. Objetivou-se espacializar os fatores físicos, climáticos e socioeconômicos que melhor discriminam a produção de caprinos leiteiros no Brasil. Foram analisados índice de produção de leite por cabra; índice de produção de caprinos; índice de produção de leite, amplitude da temperatura; temperatura média; precipitação; índice normalizado de diferença vegetativa; umidade relativa do ar; altitude; estabelecimentos agropecuários; estabelecimentos com pastagem nativa; estabelecimentos com pastagens de boa qualidade; estabelecimentos com recursos hídricos; estabelecimentos que recebem orientação técnica; estabelecimentos de agricultura familiar; estabelecimentos de agricultura não familiar e índice de desenvolvimento humano. Foram realizadas análises multivariadas para espacializar as variáveis climáticas, físicas e socioeconômicas e, assim, discriminar os Estados e Regiões brasileiras. As maiores produções de caprinos e de leite foram observadas na região Nordeste. A região Sudeste apresentou segunda maior produção de leite, seguido pelo Sul, Centro-Oeste e Norte. As médias para produtividade mostraram que as regiões Centro-Oeste e Sudeste apresentaram animais mais especializados a produção de leite. As análises multivariadas evidenciaram distinções entre clusters das regiões político-administrativas do Brasil. As variáveis climáticas foram as mais importantes para discriminar entre as regiões brasileiras. A heterogeneidade dos componentes climáticos, físicos e socioeconômicos evidenciou peculiaridades em cada região. Portanto, é preciso implementar programas de melhoramento genético animal que atendam as necessidades de cada região.

Analysis of longitudinal data of beef cattle raised on pasture from northern Brazil using nonlinear models.

This study was undertaken to aim of estimating the genetic parameters and trends for asymptotic weight (A) and maturity rate (k) of Nellore cattle from northern Brazil. The data set was made available by the Brazilian Association of Zebu Breeders and collected between the years of 1997 and 2007. The Von Bertalanffy, Brody, Gompertz, and logistic nonlinear models were fitted by the Gauss-Newton method to weight-age data of 45,895 animals collected quarterly of the birth to 750 days old. The curve parameters were analyzed using the procedures GLM and CORR. The estimation of (co)variance components and genetic parameters was obtained using the MTDFREML software. The estimated heritability coefficients were 0.21 ± 0.013 and 0.25 ± 0.014 for asymptotic weight and maturity rate, respectively. This indicates that selection for any trait shall results in genetic progress in the herd. The genetic correlation between A and k was negative (-0.57 ± 0.03) and indicated that animals selected for high maturity rate shall result in low asymptotic weight. The Von Bertalanffy function is adequate to establish the mean growth patterns and to predict the adult weight of Nellore cattle. This model is more accurate in predicting the birth weight of these animals and has better overall fit. The prediction of adult weight using nonlinear functions can be accurate when growth curve parameters and their (co)variance components are estimated jointly. The model used in this study can be applied to the prediction of mature weight in herds where a portion of the animals are culled before they reach the adult age.